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We consider four nontrivial ensembles involving Gaussian Wigner and Wishart matrices. These 
are relevant to problems ranging from multiantenna communication to random supergravity. We 
derive the matrix probability density, as well as the eigenvalue densities for these ensembles. In all 
cases the joint eigenvalue density exhibits a biorthogonal structure. A determinantal representation, 
based on a generalization of Andreief’s integration formula, is used to compactly express the r-point 
correlation function of eigenvalues. This representation circumvents the complications encountered 
in the usual approaches, and the answer is obtained immediately by examining the joint density of 
eigenvalues. We validate our analytical results using Monte Carlo simulations. 
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I. INTRODUCTION 

Wigner and Wishart matrices have been the corner¬ 
stones of random matrix theory. They find numerous 
applications in varied fields of knowledge M- The in¬ 
ception of Wigner matrices was due to Wigner who in¬ 
vestigated some special large dimensional random matri¬ 
ces to predict properties of the eigenfunctions and eigen¬ 
values of complicated quantum mechanical systems, in 
particular heavy nuclei HI- It turns out that certain 
spectral characteristics of these matrices, such as semicir¬ 
cular distribution of eigenvalues, are universal and in fact 
shared by a wider class of matrices which are now collec¬ 
tively referred to as Wigner matrices. See, for example, 
Refs. @, 0 for recent surveys. An important family of 
Wigner matrices is realized when the matrix elements are 
taken as Gaussian random variables. The resulting en¬ 
sembles are referred to as Gaussian ensembles 0,01- In 
the present work we use the terms Wigner and Gaussian 
Wigner interchangeably to mean complex Wigner matri¬ 
ces with Gaussian entries, more specifically matrices from 
the Gaussian unitary ensemble |l, 1, 0- Wishart matri¬ 
ces predate even Wigner matrices and have their origin 
in the field of multivariate statistics. These were intro¬ 
duced by Wishart who derived the generalized product- 
moment distribution in normal multivariate population 
samples 0]. This distribution is now referred to as the 
Wishart distribution. In what follows, we will be con¬ 
cerned with ensembles comprising complex Wishart ma¬ 
trices. 

While Wigner and Wishart matrices themselves offer 
plenty of aspects to explore, interestingly, various com¬ 
binations of these also turn out to be of crucial impor¬ 
tance. Many such matrix ensembles have their origin in 
the area of multivariate statistics jlol - [0 ]. A classic ex¬ 
ample is the Jacobi (MANOVA) ensemble which incorpo¬ 
rates two Wishart matrices in a nontrivial manner, and 
arises in the problems of quantum conductance [13Hl8j 
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and multiple channel fiber optics communication 0,0. 
Remarkably, it also pops up in something as remote as a 
microscopic model of bus transport system 21]. In recent 
years several other important matrix models have been 
explored. Some notable examples include product of 
complex Ginibre matrices [22-28[, Cauchy-Lorentz 
|33| , sum involving Wigner and Wishart matrices [U 
and product of truncated unitary matrices 0] . In addi¬ 
tion to their natural connection with multivariate statis¬ 
tics, these are of interest to the fields of telecommunica¬ 



tion [22t l26ll38ll. finance [41|, |42|, and random supergrav¬ 


ity theory [3443 


In the present work we proceed further in exploring 
such exotic ensembles and consider four important ma¬ 
trix models involving Wigner and Wishart matrices. The 
first one comprises a ratio involving two Wishart matri¬ 
ces, the second one consists of the weighted sum of a 
Wigner matrix and a Wishart matrix, the third is the 
product of a Wigner matrix and a Wishart matrix, and 
the fourth one embodies the weighted sum of two Wishart 
matrices. We derive the probability density function for 
these matrices, and then work out the eigenvalue statis¬ 
tics. The joint density of eigenvalues for these matrix 
models exhibit biorthogonal structure. A determinantal 
representation, based on a generalization of Andreief’s 
integration formula 0 - 0 , is used to compactly express 
the r-point correlation function for all these ensembles. 


The rest of the paper is organized as follows. We 
start with a brief discussion of biorthogonal ensembles 
in Sec. im and present the result for r-level correlation 
function for eigenvalues. Sections [TTH EJ are devoted to 
the exact results for the above mentioned matrix ensem¬ 
bles which involve nontrivial combinations of Wigner and 
Wishart matrices. We conclude in Sec. ivTn with a brief 
summary and outlook. Some relevant derivations are pre¬ 
sented in the Appendices. 
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II. BIORTHOGONAL ENSEMBLES 

Biorthogonal ensembles arise naturally in the study of 
eigenvalue statistics of two matrix models [Ifl [47j . More¬ 
over, matrix ensembles with a unitary invariance break¬ 
ing external source also give rise to biorthogonal struc¬ 
ture gmg. These ensembles exhibit rich mathemat¬ 
ical structure and, at the same time, find applications 
in several important problems which range from quan¬ 
tum transport to multiple antenna telecommunication, 
to two-dimensional gravity [2214251 15 D - I 57 ] ]. 

We are interested here in biorthogonal ensembles of the 
Borodin type [HI, which possess the following structure 
for joint density of its eigenvalues {A} (= {Ai,..., An}): 


correlation function in terms of a determinant of a r- 
dimensional matrix with entries containing certain two- 
point kernel. However, it requires inversion of a matrix. 

In the following we use a generalization of Andreief’s 
integration formula to express the r-point correlation 
function in terms of the determinant of an (n + r)- 
dimensional matrix [HI, |45[ : 

r 

i2r(A 1 ,...,A r ) = (-l) r n!C , n^(A,) 

1=1 

x k=l,...,n fa) 

[fj(Xk)]j-l,...,n [hj,k]j=l,...,n 

k=l,...,r k=l,...,n 


P({A}) = C'A n ({A})n«'(A«) • ./)bV!} ; .,-,.„ ■ (1) 

1=1 

Here w{ A) is a well-behaved weight function in the desired 
domain, and |_| represents determinant. Also, A„({A}) 
is the Vandermonde determinant, 

A„({A}) = |Af 1 |,-, fe= i,...,„ = n(A, - A fc ). (2) 

j>k 

The normalization factor C follows by expanding the de¬ 
terminants and performing the integrals. The ensuing 
expression can again be represented as a determinant, as 
asserted by Andreief identity [43[ . We have 

C = Tl\ \hj,k\j,k=l,...,n: (3) 

where 

hj,k = J dXw(X)fj(X) A fe_1 . (4) 

For the special case of fj(Xk) = A(T , we obtain the joint 
probability density of eigenvalues for a unitary random 
matrix ensemble. We note that if we replace A„({A}) 
by some other determinant \gj( Afc)|, then we have the 
most general form of biorthogonal ensemble, as defined 
by Borodin (54| • The approach for calculating correlation 
function, as discussed below, extends to these as well. 

We would like to remark that the biorthogonal ensem¬ 
ble of Borodin type emerges after integrating out one 
set of eigenvalues (corresponding to one of the matrices) 
from the joint probability density of eigenvalues for two- 
matrix model; see for example Appendix [Bj 

The r-point correlation function (1 < r < n) corre¬ 
sponding to Eq. Q is defined as [I| 

Rr( Ai, A r ) = ———— J dX r+ i J -P({A}). 

(5) 

The evaluation of this correlation function usually relies 
on the explicit construction of biorthogonal polynomials. 
In [54} a recipe has been provided to write down the 


In the above expression 0 represents r x r block with all 
entries 0, and fj(Xk), hj t k are as appearing in Eqs. <[T|) 
and CD , respectively. In Appendix [A] we provide a proof 
of Eq. © based on mathematical induction. The above 
representation for correlation function altogether circum¬ 
vents the complications encountered in the approaches 
described above, and an explicit answer is obtained at 
once. For small n, r Eq. 0 is advantageous in the sense 
that it can be readily implemented and evaluated in com¬ 
putational packages such as Mathematica [58}. In par¬ 
ticular the first-order marginal density of eigenvalues, 
which is related to the one-point correlation function as 
p{ A) = I?i(A)/n, is given by 


pW 


~(n 


l)!Cw(X) 


0 [A fc - 1 ] fe= i.„ 

[fj (A)]j=l,...,n [hj,k]j= l,...,n 


A similar form has been used in [37i . Ivil to express 
the marginal density of eigenvalues. On the other ex¬ 
treme, if we consider r = n, then the determinant in 
Eq. © collapses to the product of determinants |A^ _1 | 
and |/j(Afc)|, along with the factor (—1)", and thereby 
yields n!P({A}), as expected. 

As discussed in the Introduction, in the following 
sections we consider four matrix ensembles where such 
biorthogonal structure emerges. The joint density of 
eigenvalues for these ensembles appear in the form of 
Eq. ©, and hence the r-point correlation function can 
be written down immediately with the aid of Eq. ©. 


III. RATIO INVOLVING TWO WISHARTS 

A. Matrix model and probability density 

We consider an ensemble oinxn dimensional complex 
matrices 


H = {aA)(t n + bB)-\ (8) 

where a and b are some non-negative scalars (for defi¬ 
niteness), and A and B are positive-definite Hermitian 
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matrices, respectively, from the complex Wisliart distri¬ 
butions 

Va{A) oc e- trA \A\ nA - n , V b (B) oc e~ tr B \B\ nB ~ n . (9) 

Here ua,ti b > n are the respective degrees of freedom 
for the two distributions. We may refer to the ensem¬ 
ble described by Eq. © as a quotient ensemble. For 
b —> 0 we have the usual complex Wishart, while the limit 
a = b —> oo leads to the ensemble AB -1 , which is a multi¬ 
variate generalization of the beta distribution of the sec¬ 
ond kind [60}. We also note that (l ra +bB)~ 1 / 2 (ciA)(t n + 
bB )~ 1//2 , aA(t n + bB) -1 and (1 n +bB)~ 1 (aA) share the 
identical nonnegative eigenvalues as they correspond to 
the same generalized eigenvalue problem and lead to the 
secular equation \aA — X(t n + bB)\ = 0. We will see be¬ 
low that the above construction leads to a very interest¬ 
ing matrix model whose probability density involves con¬ 
fluent hypergeometric function of the second kind (Tri- 
comi’s function) with matrix argument a Moreover, 
this matrix model is of direct relevance to the problem 
of multiple antenna relay systems [62J. 

The probability density of H can be calculated as 

V h (H) = J d[A]P A (A) I d[B]V B {B) 

x6(H - (aA)(t n + bB)~ 1 ). (10) 

Here the delta function with matrix argument represents 
the product of delta functions with scalar arguments, one 
for each independent real and imaginary component of 
H — (aA)(l n + bB) -1 . Also, d[A], etc. represent the 
flat measure involving the product of the differentials of 
all independent variables occurring within the matrix. 
Implementation of the Fourier representation for delta 
function and the cyclic invariance property of trace gives 

V h (H) cx J d[K\ J d[A]P A (A) j d[B\P B (B) 

xe itr KH e ~iti (aAtln+bB )- 1 K) 


The matrix K in the above equation possesses symmetry 
properties identical to those of H — (aA)(t n + bB) -1 . 
Using Eq. ©, reordering the integrals, and considering 
the transformation K -A- ( l n + bB)K, we obtain 

V h {H) oc J d[B]e~ trB \B\ nB ~ n \t n + bB\ n 

x J d[K}e itrKH{1 "- +bB) J d[A}e- tlA ^ +iaK) \A\ nA ~ n . 

( 12 ) 

Integration over A yields 

V h {H) oc J d[B]e~ trB \B\ nB ~ n \t n + bB\ n 

x [ d[K]e itlKH( ' ln+bB) \t n +iaK\~ nA . (13) 


The K integral can be identified as a variant of Ingham- 
Siegel-Fyodorov integral [63j] and leads to 

V h {H) oc J d[B}e- tlB \B\ nB ~ n \l n + bB\ n 

xe -tra- 1 H(t n +bB)^ H ^n A -n^ tn + n A -n _ 

We may write 

V h {H) oc e- a ~ ltlH \H\ nA ~ n A>(H), (15) 

where 

$(Ff) = f d[B}e~ tT( ‘ ln+a ~ lbH ' >B \B\ nB ~ n \t n + bB\ nA . 

(16) 

&(H) can be expressed in terms of the confluent hyperge¬ 
ometric function of the second kind (Triconri’s function) 
with matrix argument [6lJ, 

vk(a, 7 ;X) = -^ f d[Y]e- tlXY \Y\ a ~ n \l n +Y\'y- a - n , 
Tn(a )J 

(17) 
as 

= frlne + UB + n\ (b~ l l n + a -1 H)). 

(18 ) 

Here T n (riB ) is the multivariate Gamma function: 

n 

r n (a)=n^-^l[T(a-j + l). (19) 

l=i 

Thus, we finally have the result 

V h (H) oc e - a ~ ltrH \H\ n A - n 

x >F(riB, riA + riB + n; (b -1 t n + a -1 H)). (20) 


B. Eigenvalue statistics 


We now derive the joint density of eigenvalues for the 
matrix model of Eq. ©. As implied by the result in Ap¬ 
pendix!© ik(a, 7 ; X) of Eq. (fT71) admits the following de- 
terminantal representation in terms of elements involving 
hypergeometric function of the second kind (Tricomi’s 
function) with scalar argument [64l 1651]: 


ik(a, 7 ; X) oc 


1 


A((4) 


\U{a-j+l, 7 -j-n+ 2 ; a;*:) |j,fc=i,.. 


( 21 ) 

Here ccy’s are the eigenvalues of X. The joint density of 
eigenvalues (0 < Ai,A„ < oo) of H , therefore, follows 
immediately from Eq. m, and possesses the biorthogo- 
nal structure as in Eq. © with 


w{ A) = e- x / a X nA ~ n , (22) 


fj( A*) = C7( 


n B ~ j + l,n A + n B 
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FIG. 1: Eigenvalue densities for the ratio of Wisliart matrices, Eq. (f8|) . (a) Two point correlation function for 
n = 3, ua = 20, ns = 21, a = 2, b = 1/5; (b) marginal density for n = 4, ha = 14, ns = 9, a = 1 and different b 
values, as indicated. The symbols (circles, squares, triangles) shown in (b) are using numerical simulation. 



FIG. 2: Marginal densities for two Wishart matrices 
(Wishart A , Wishart B) and the ratio, as defined in 
Eq. ©. The parameters are n = 6, ua = 9, ns = 18, 
and a = 2, b = 1/2. 


The hj^k of Eq. (HI) is obtained as 


not shown here for the sake of clarity, a two-dimensional 
histogram obtained using Monte-Carlo simulation agrees 
well with the analytical plot. In Fig. marginal density 
of eigenvalues is shown for parameter values mentioned 
in the caption. In this case simulation results are also 
depicted with the aid of symbols, and are in excellent 
agreement with the analytical curves. As already indi¬ 
cated, the parameters a and b give freedom to realize a 
variety of densities using two Wishart matrices, the exact 
outcome being dependent on the dimensions of the con¬ 
stituent matrices. As an example, in Fig. [2] we show the 
density corresponding to the quotient ensemble defined 
by Eq. © along with the densities of the constituent 
Wishart matrices, which can be calculated using the re¬ 
sult 


Pw(A) = 


F(n) 


-A s e“ A 


T(s + n) 

x [L^A^+^A) - L«(A) J L,t + i 1) (A)]. (26) 


h jjk = a nA ~ n+k T(n A - n + k) 

x u(n B — j + l,n B + n - j - k + 2; (24) 

We used here the integral result 

OO 

dz z c e~ z U(a,b; z + m) = T(c+ l)U(a,b — c — l;m), 

(25) 

which holds whenever the integral is convergent. There¬ 
fore, r-point correlation function and the marginal den¬ 
sity follow immediately from Eqs. © and 0. 

In Fig. [la| we show the two-point correlation function 
for parameter values indicated in the caption. Although 


Here L^\ A) represents associated Laguerre polynomial 
of degree fi, and the parameter s is given by ua — n 
or ub —n, i.e., it is the difference of degree of freedom 
and dimension of the Wishart matrix. We should un¬ 
derline that the resulting spectra, tunable by a and b , 
can have a crucial role in deciding the behavior of met¬ 
ric which follow from the eigenvalue statistics, such as 
channel capacity and outage probability in the case of 
multiple access channel (MAC) and interference channel 
(IC) in multiple-input multiple-output (MIMO) commu¬ 
nication jfi'j . 
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FIG. 3: Eigenvalue densities for weighted sum of Wigner and Wishart matrices, Eq. m- (a) Two point correlation 
function for n = 3, ns = 4, a = 4, b = 1; (b) marginal density for n = 5, Ub = 7, b = 1 — a and a = 9/10,6/10, 3/10 as 
indicated in the figure. 



FIG. 4: Marginal densities for Wigner matrix, Wishart 
matrix and their weighted sum, as given in Eq. m- 
The parameter values are n = 6, ns = 9, a = 1, b = 1/2. 


IV. WEIGHTED SUM OF A WIGNER AND A 
WISHART 


sian unitary) ensemble. On the other hand, for a —> 0, 
with b > 0, we obtain the Wishart (Laguerre unitary) 
ensemble. Therefore, by considering b = 1 — a, and by 
varying a between 0 and 1, we have an ensemble which 
interpolates between the Wishart and Wigner ensembles. 
To the best of our knowledge, for this matrix model only 
the first order marginal density of eigenvalues is known 
in the large n asymptotic regime using the tools of free 
probability f66j . A matrix ensemble similar to that in 
Eq. (l27l) has been used to model the Hessian matrix in 
the context of supergravity (34| - [36j . 

To obtain the probability density function for H we 
introduce the Fourier representation of delta function as 
in Eq. (HD. Reordering of the integrals, and use of the 
cyclic invariance property of trace then gives 


V h (H) a 
x 




i tr (H—bB)K 


(29) 


Evaluation of the Gaussian integral involving A leads to 


A. Matrix model and probability density 

We now consider an ensemble comprising weighted sum 
of Wigner and Wishart matrices: 

H = aA + bB. (27) 

Here A and B are respectively n-dimensional Hermitian 
and positive-definite-Hermitian matrices from the distri¬ 
butions 

Va(A) (xe~ trA \ V b {B) oc e~ trB \B\ nB ~ n , (28) 

and a, b, as before, are non-negative scalars. Also, ub > 
n. For b —> 0, with a > 0, we have the Wigner (Gaus¬ 


P h (H) oc J d[B)e- trB \B\ nB - n 

x J d[K}e-£ tlK2 e itr(H - bB)K . (30) 

The Gaussian integral over K can also be performed and 
yields 

V h (H) oc e"^ trff2 4>(tf), (31) 

where 

$(F0 = J d[B}e~ tr b 2 e tT ^ H ~^ ln ^ B \B\ nB ~ n . (32) 
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B. Eigenvalue statistics 

We now calculate the eigenvalue statistics correspond¬ 
ing to Eq. m- Using the result in Appendix [B] we 
know that Q(H) is determined solely by the eigenvalues 
(—00 < Ai, X n < oo) of H as 

Q(H) OC (33) 

I 


where 

fj{X h ) = [°° (34) 

Jo 

This integral can be evaluated in terms of confluent hy¬ 
pergeometric function of the first kind (Rummer’s func¬ 
tion) [(34l and leads to the joint density, Eq. (JT]) , 
with m 


fA Afc) = 


n B ~ j + 1 


1 E 1 


n B ~ j + 1 1 / A fc 


+ 


A k 

a 


a 

2b 


n B — j + 2 


i F 1 


2 ’ \ a 
n B — j + 2 3 


A k 

a 


a 

2b 


(35) 


The weight function is read from Eq. m as 


W (A) =e- A2 /“ 2 . (36) 

In this case obtaining a closed form for hj k requires some effort. A possible representation is in terms of hypergeometric 
2 F 2 Hill: 


x/Trb nB ~ j+k 

hj,k = -r—■-r(n B - j + k ) 2 F 2 


1 — k 2 — k 1 — n B + j — k 2 — n B + j — k a 2 

2 ’ 2 ’ 2 ’ 2 ’IP 


(37) 


With the above explicit results, the r-point correlation function of Eq. © is readily obtained. 


We show the two-point correlation function surface in 
Fig. I3al The marginal density is shown along with the 
Monte-Carlo simulation outcome in Fig. I3bl In particu¬ 
lar, for Fig. [3b] we have considered b = 1 — a. Therefore, 
a crossover is seen from Wigner density (semicircle type) 
to Wishart density (Marcenko-Pastur type). In Fig. Q] 
we compare the eigenvalue density of the composite en¬ 
semble with the eigenvalue density for the constituent 
Wishart ensemble given by Eq. (l26l) and that of the Gaus¬ 
sian Wigner ensemble evaluated using 


— X 2 

Pgw(A) = e [g n (A)g n (A) - H n _i(X)H n+ i(X)]. 
2 n v '7rn! 

(38) 


Here A) represents the Hermite polynomial of degree 
p,. We can see that the Wishart constituent of the com¬ 
posite matrix tends to keep the eigenvalues in the positive 
half of the real line, while the Wigner part pulls them to¬ 
ward the negative half and tries to make the density sym¬ 
metric about zero, thereby giving rise to an interesting 
hybrid density. 


V. PRODUCT OF A WIGNER AND A 
WISHART 

A. Matrix model and probability density 

We now consider an ensemble defined by 

H = AB, (39) 

where A and B , respectively, are Wigner and Wishart 
matrices from the distributions given in (1281) . We note 
here that A is Hermitian but B is Hermitian and positive- 
definite as well; therefore, the signs of eigenvalues of 
H are decided by the respective signs of eigenvalues of 
A [67]. As a consequence we expect the resultant first or¬ 
der marginal density to be symmetric about the origin, 
similar to that in the Wigner case. 

We introduce the matrix delta function, as in Eq. CD. 
to obtain 

V h (H) oc J d[K] j d[A] J d[B}e itr K ^ H ~ AB ) 

x e- trA2 e- tlB \B\ nB - n . (40) 

We reorder the integrals and use the cyclic invariance of 
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FIG. 5: Eigenvalue densities for product of Wigner and Wishart matrices, Eq. (l39l) . (a) Two point correlation 
function for n = 3 ,ns = 5; (b) marginal density of eigenvalues for n = Hb = 2 and n = 7,ns = 8. 



FIG. 6: Marginal densities for Wigner matrix, Wishart 
matrix, and their product as in Eq. (1391) . The dimen¬ 
sions of the matrices are given by n = 6, ns = 9. 


The A'-integral can now be performed |j33[ and yields 

V h (H) oc \H\ nB - n §{H) (44) 

with 

$(i?) = J d[^]e“ tr(A2+ ^ lff) |A|- riB QiA^H). (45) 

Here 0(_) represents the Heaviside theta function and 
requires A~ 1 H to be positive-definite for a non-vanishing 
result. 


B. Eigenvalue statistics 

With a little modification the result in Appendix iBlim- 
plies that 4>(Ff) is determined by the eigenvalues (—oo < 
Ai,..., A n < oo) of H as 


trace to get 

V H {H)ot J d[A}e~ tlA2 J d[K]e itrKH 

x J d[B]e- tlB( -^ +iKA) \B\ nB - n . (41) 

Integral over B can be done to give 

Th{H) oc J d[A]e~ tlA2 J d[K]e itlKH \t n + iKA\~ nB . 

(42) 

We now employ the transformation K —> KA 1 , which 
leads to 

V h (H) oc J d[A]e~ tr a 2 \A\~ n 

[ d[K]e itlKA ~ lH \l n + iK\- nB . (43) 


m) 


A({A}) 




(46) 


where 


fj{ A fc )= / dfif,- nB+n+j - 2 e-^~ Xk ^, (47) 


with u = —oo for A < 0 and u = oo for A > 0. This 
integral can be evaluated compactly in terms of Meijer 
G-function 64] as 


—UB+n+j — 1 


/.*(**) = 77f 


2 -ns+n+j 


X Gq;° UB - n - j + 1 ns ~ n - j + 2 




X 


(48) 

























Meijer G-functions have also appeared in the correlation 
kernels for product of complex Ginibre matric es |22 h 28[ , 
and product of truncated unitary matrices [401] . The 
weight function w(X), in view of Eq. fill) . is 


w(X) = X nB ~ n , (49) 

which leads to the following expression of h.j.k'- 

hj,k = 1 + ( : lp+fc r (n B -n + k) t( • (50) 


been considered in {38| and the results used for the analy¬ 
sis of multiuser communication employing multiantenna 
elements such as multiple-input multiple-output (MIMO) 
multiple access channel (MAC). 

For the matrix model of Eq. ED, with parameter m = 
riA + n,B — n, the probability density function satisfied by 
matrix H reads [37| 

V h {H) oc |iZ| m e -tr ( a_1 ‘ f/ ) 

x 1 F 1 (n B -n A + n B -(a^ 1 t n -b~ 1 T,~ 1 )H), (53) 


Equation © now determines correlation functions of all 
orders for the matrix model ED- 

Figure [Sa] shows the two-point correlation function of 
eigenvalues, while Fig. [5b] depicts the marginal density. 
For n = n B the density exhibits a logarithmic singular¬ 
ity at A = 0. This can be seen in n = n B = 2 plot 
in Fig. [5b] In Fig. [ 6 ] we display the eigenvalue density 
for the product of Wigner and Wishart matrices, along 
with the densities for the constituent matrices calculated 
using Eqs. (1551) and (1551) . Figure [ 6 ] should be compared 
with Fig. [4] where we have used matrices with dimensions 
same as in the present case. The distinct nature of the 
resultant densities in these two cases is expected because 
of very different characteristics of the underlying compo¬ 
sition. The shapes of the marginal density curves here 
are reminiscent of the density of eig envalues of adjacency 
matrices in scale free networks [68H7CM and matrices de¬ 
fined on Poissonian random graphs [7l| . 


where \ F\ is confluent hypergeometric function of the 
first kind (Rummer’s function) with matrix argument: 


iFi(a,7;X) = 


1 


B n (a, 7 — a) 


x / d[Y}e tlXY \t n -Y\ a - n \Y^- a - n . 

Jo 

Here B n (a, 7 ) is the multivariate beta function: 


(54) 


B n (a, 7 )= / d[Y}\t n -Y\ a - n \Yy- n . (55) 
Jo 


Similar to the beta function with scalar arguments, it is 
related to multivariate gamma function in Eq. ED as 


B n (a, 7 ) 


F n (q) r n ( 7 ) 
r„(a + 7 ) 


(56) 


VI. WEIGHTED SUM OF TWO WISHARTS 


B. Eigenvalue statistics 


A. Matrix model and probability density 


The joint probability density of eigenvalues (0 < 
Ai,..., A„ < oo) for Eq. ED is given by Eq. © with 


We finally consider the matrix model 


w( A) = X m e~^ a , (57) 


H = aA + bB , (51) 

where A and B are n-dimensional positive-defmite- 
Hermitian matrices, respectively, from the distributions 

P{A) oc e _tr A |A| raA_n , p(B) oc e ~ tli: ~ lB \B \ nB - n , 

(52) 

with n a , n B > n. The parameters a, b are again non¬ 
negative scalars. We have taken the covariance matrix 
equal to identity matrix for A, while for B we have 
assumed an arbitrary (positive definite) covariance ma¬ 
trix. The latter constitutes the correlated variant of the 
Wishart ensemble. The above matrix model has been 
considered in [37| and exact results have been obtained 
for the matrix probability density, the joint probability 
density of eigenvalues, as well as the first order marginal 
density. 

We would like to remark that if one considers covari¬ 
ance matrices proportional to identity matrix only, then 
the problem can be solved for the weighted sum of arbi¬ 
trary number of Wishart matrices. Such a scenario has 


fj{ Afc) = iF 1 (n B - n+ 1; m + 1; - ^-)a(58) 

where <jj are the eigenvalues of E s 37j|. Also, hjj- can be 
obtained using the result 


J^dX X m e~ sX 1 F 1 (a ; b; cX) = 2*1 (a; m+ 1; 6 ; J), 


valid for convergent scenarios, as 
h j>k = a m+k r (to + k) 
x 2 A] ( n B — n+l;m — n + k; m + l; 1 ~ 


ba 3 


(59) 


■ (60) 


Consequently, we obtain an explicit result for the r-point 
correlation function. 

Figure [7a] shows the two-point correlation function of 
eigenvalues for matrix model given in Eq. ED, while 
Fig. [7b] depicts the marginal density. In Fig. [8] we 
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FIG. 7: Eigenvalue densities for weighted sum of Wishart and correlated-Wisliart matrices, Eq. m- (a) Two point 
correlation function for 6=1; (b) marginal density for 6=1/3 and 1. Common parameters for both the figures are 
a = 1/4, n = 4, ha = 10, ns = 11, and (cri, cr 2 , ct 3 , 0 - 4 ) = (5/2,1/3, 2, 7/4). 



FIG. 8: Marginal densities for Wishart matrix, 
correlated-Wishart matrix, and their weighted sum as 
given in Eq. ED- The parameter values are 
n = 6, tia = tib = 9 , a = 6 = 1 , and 
( CTl ,..., <76) = (4,20/3,5/2, 11 / 9 , 4/3, 7/8). 


compare the densities of Wishart, correlated-Wishart 
and their weighted sum. For density of the correlated- 
Wishart we have used the following result, which also 
follows with the aid of Eq. 0: 


Pew (A) = — 
n 



X s 


A(w)n;u^ +lr u+s) 

1 ,...,n [r(fc + s)o| +s ]j = i,... iri 

fe=l,...,n 


(61) 


where s = ns — n. We note that while matrix models (|8|) 
and m can lead to similar densities in particular sce¬ 
narios (e.g. when 6 —> 0), in general they exhibit different 
behavior. For instance, the present model can not lead to 


the nmltivariate-beta-distribution kind of density which 
can be achieved using Eq. 0 for a = b —> 00 . 


VII. CONCLUSION AND OUTLOOK 


We considered four important matrix models which 
lead to biorthogonal structure in their joint eigenvalue 
densities. These matrix ensembles play important roles 
in several areas, which range from multiple antenna com¬ 
munication theory to supergravity theory. We evaluated 
the matrix distribution, as well as the joint eigenvalue 
density for these ensembles. With the information of 
joint density, we also presented determinantal expres¬ 
sion for eigenvalue correlation function of arbitrary or¬ 
der. This representation follows from a generalization 
of Andreief’s integration formula. Since knowledge of 
the correlation function gives access to the prediction of 
statistical behavior of observables of interest in a given 
problem, we believe that the exact results derived here 
will find interesting applications in several fields. 

As continuation of the present work, an immediate di¬ 
rection to pursue could be the investigation of the behav¬ 
ior of extreme eigenvalues of the composite matrix mod¬ 
els, and its comparison with that of the extreme eigen¬ 
values of the constituent matrices. This will give a better 
insight into the mechanism by which the redistribution of 
eigenvalues takes place. Since all the matrix models con¬ 
sidered here possess biorthogonal structure in their joint 
eigenvalue density expressions, exact results are possible 
for the gap probabilities and densities of extreme eigen¬ 
values |60|. 

While we considered here ensembles comprising com¬ 
plex matrices, the cases of real and quaternion matrices 
are also important and can be explored. However, solv¬ 
ing these ensembles poses serious challenges because of 
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unavailability of group integral results similar to those 
in the case of unitary group. Another interesting direc¬ 
tion can be the analysis of the spectra of the composite 
matrices in large dimension limit, and to look for univer¬ 
salities. 


Appendix A: Correlation function 

We will use mathematical induction to prove Eq. (1(11) . 
Equations © and © are defined for r = 1,2 |75| . 
From the definition of correlation function, Eq. ©, it is 
clear that 

R r —i (Ai,..., A r _i) ■ “ / dX r R r (Ai, • ■ •, A r ). 

n — r + 1 J 

(Al) 

For r = n Eq. © clearly holds, since in this case the 
determinant in Eq. © factorizes into the product of two 
determinants and produces n!P({A}). Let us assume it 
is correct for r = s. We will prove that given this, Eq. © 
holds for r = s — 1 as well. 


Using Eq. (ED we obtain 


Ps-i(Ai,..., A s _i) = 


(—1 ) s n\C 
n — s + 1 




i=i 


We expand the determinant using the s’th row: 


r 

IP]j=i,...,s i^j 1 ]j=i>—.* 

/ dX s w(X s ) 

k=l,...,s k=l,...,n 

[/j(Afc)]j=l,...,n [h'3,k\j— 

k=l,...,s k=l,...,n 


(A2) 


[0] J=1,...,S 

A') ., 

n 

[0]j=l,...,s-l [rf 1 ]j=l,...,s-l 

k=l,...,s k=l, ...,n 

k=l,...,s k=l,...,n 

[fj (Afe)]j=l,...,n [hj,k]j= 

= Et-i^Ar 1 

[/j'(Afe)]j'=l,...,n [hj.k\j=l,...,n 

k=l,...,s 

k=l,...,n 

M=1 

k=l,...,s k=l,...,n 

(kjtfi) 


(A3) 


We now insert the w( A s )A(f 1 in the s’th column, and perform the A s integral. Using the definition of h h k given in 
Eq. © , we obtain 


XU 1 ) 1 

fl=l 



,o, 

to. 

II 

In 

1 

a'; 'j i., i 



k=l,...,n 

[/j'(Afc)] j=l,...,n 


[hj,k]j=l,...,n 

k=l,... ,s —1 


k=l,...,n 

(kjtii) 


(A4) 


Performing separate row interchanges in the determinants appearing in the sum, we arrive at 


n 

[0]f=l,...,s-l 

[(1 - <W)A ) 1 ]i=i,..., s -i 

11=1 

k=l,...,s —1 

[/j(Afc)] j=l,...,n 
k=l,...,s—l 

fc= 

[hj,k]j=l,...,n 


where S^ is the Kronecker-delta function. 


(A5) 


Using multilinearity property in first s — 1 rows in de¬ 
terminant appearing in each of the terms in the above 
summation, we find that it gives rise to 


(-1) 1 (n-s + 1) 


[0]j=i.. 

k=l,...,s-l 

[fj(^k)] j=l,...,n 
k=l,...,s-l 


vfc-li 


[A‘- 

fc=l,...,n 

[hj,k\j=l,...,n 




(A6) 


Plugging this back in Eq. (IA2I) . we obtain an expression 
for i? s -i(Ai,..., A s _i) which is consistent with Eq. ©, 
and hence the desired result follows. 


Appendix B: Matrix Integral 


Consider n-dimensional Hermitian matrices X and Y. 
We are interested in evaluating integral of the form 


$(X) = 


d[Y } e~ stlXY F{Y), 


(Bl) 


where s is a scalar and F(Y) is a unitarily invariant ex¬ 
pression involving Y, such that the above integral is con¬ 
vergent. We note that Eq. EH) is a matrix generalization 
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of Laplace transform. If x and y are the diagonal matri¬ 
ces consisting of eigenvalues of X and Y, then 


$(X) = / dy-i ■■■ dy n A 2 n ({y})F(y ) 

Jo Jo 

x [ dn{U)e~ sU{ - xU] yU) , (B2) 

JUn 


where dy(U) represents the Haar measure over the group 
lA n of ?r-dimensional unitary matrices. The unitary group 
integral can be performed using the celebrated Harish- 


Chandra-Itzykson-Zuber formula [72;, 173] and leads to 

oc A (t \\ / ri 2h ■ dy n A n ({y}) 
A„({a;}) J o J o 

xF(y)\e- s ^\. k= i_ n - (B3) 

Now if F(y) is expressible in terms of certain weight func¬ 
tions u(yj) as F(y) = nj=i u {Uj)i then integral over y 
can be performed and results in 

. (B4) 

where 

/*00 

fj(xk)= / dyu{y)y J ~ 1 e~ SXkV . (B5) 

Jo 

Note that we may consider j —> n — j + 1 (or/and k —> 
n — k + 1 ) for fj{xk ) within the determinant in ([I]) and 
then accordingly modify the rest of the results in Sec. ITT1 
which depend on fj(xk)- 
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